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Abstract 

Spectroscopic techniques are very essential tools in studying electronic structures, 
spectroscopic constants and energetic properties of diatomic molecules. These tech¬ 
niques are also required for parametrization of new method based on theoretical analy¬ 
sis and computational calculations. In this research, we apply the proper quantization 
rule in spectroscopic study of some diatomic molecules by solving the Schrodinger 
equation with two solvable quantum molecular systems-Tietz-Wei and shifted Deng- 
Fan potential models for their approximate nonrelativistic energy states via an appro¬ 
priate approximation to the centrifugal term. We show that the energy levels can be 
determined from its ground state energy. The beauty and simplicity of the method 
applied in this study is that, it can be applied to any exactly as well as approximately 
solvable models. The validity and accuracy of the method is tested with previous 
techniques via numerical computation for H 2 and CO diatomic molecules. The result 
also include energy spectrum of 5 different electronic states of NO and 2 different 
electronic state of IC1. 
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1 Introduction 


The exact solutions of solvable quantum potential models have received much interest since 
they provide us some insight into the physical problem under consideration. Over the 

^E-mail: fbjamesll@physicist.net; babatunde.falaye@fulafia.edu.ng 
•[Corresponding author 

2 E-mail: sameer.ikhdair@najah.edu; sikhdair@gmail.com. 

^E-mail: majid.hamzavi@gmail.com 


1 




past years, various eigensolution techniques have been proposed to solve quantum potential 
models. Few of these methods are: formula method [I], Nikiforov-Uvarov method [2], the 
asymptotic iteration method |3j, the supersymmetric quantum mechanics [4], the factoriza¬ 
tion method [5], wave function ansatz method [6], the generalized pseudospectral method 
(GPS) [7j and the exact quantization rule (EQR) [8l (9l HIT]. Notes on these techniques can 
be found in ref. m 

Recently, the EQR has been proposed to solve the wave equations with some exactly or 
approximately solvable quantum potentials for their energy eigenvalues and wave functions. 
Nevertheless, such solutions involve highly complicated integral calculations, in particular 
when calculating the quantum correction term. Therefore, in order to avoid these difficulties, 
Serrano et al have proposed a new way to treat these problems and called it the proper 
quantization rule (PQR) [12) • Furthermore, it has been shown that PQR is a powerful tool 
in finding the eigenvalues for all solvable quantum systems [131 HU- 

Furthermore, the study of the bound state processes is fundamental to understanding the 
molecular spectrum of the diatomic molecules and their properties in quantum mechanics. 
In light of this, there has been a growing interest in searching for the empirical potential 
functions for diatomic molecules in chemical physics and related areas [EKE]- The reason 
is that such potentials provide the compact way to summarize what we know about a 
molecule. Thus, efforts to construct a universal potential function that fit experimental 
data in computational chemistry have been made by many researchers. It has been found 
that the exponential type molecular potentials are better than the harmonic oscillator in 
simulating the atomic interaction for diatomic molecules. 

In this context, to achieve the goal of the present work, we study the spectrum of some 
diatomic molecules using two exponential-type of molecular models; namely, the Tietz-Wei 
and shifted Deng-Fan potential models PHEJ. The bound state solution of the Schrodinger 
equation with these diatomic molecular potentials provides the rotational-vibrational energy 
states of the diatomic molecules in an accurate manner. We apply PQR to obtain the energy 
spectrum of the two molecular potential models and then obtain the rotational-vibrational 
energy states for various diatomic molecules. 

This work is organized as follows. In section 2, we briefly introduce the PQR. In section 
3, we apply the method to obtain the energy spectrum of the Schrodinger equation with 
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Tietz-Wei and shifted Deng-Fan molecular potentials. We give our numerical results and 
discussions in section 4. Some concluding remarks are given in section 5. 


2 A brief review to proper quantization rule 


In this section, we give a brief review to this method 0 m H33. The one dimensional 
Schrodinger equation takes the form: 


/A 2 ii 

— ^(a:) + ~^l E Mx) = 0, 


( 1 ) 


and can be re-written as 


<f>'{x) + 4>(x) 2 + k(x) 2 = 0, with k(x) = ^j^[E — V(x)\, 


( 2 ) 


where (j){x) = ij)'(x)/^(x) is the logarithmic derivative of the wave function x ). The 
prime denotes the derivative with respect to the variable x. n denotes the reduced mass 
of the two interacting particles. k(x) is the momentum and V(x) is a piecewise continuous 
real potential function of x. According to Yang [IS] “For the Sturm-Liouville problem, 
the fundamental trick is the definition of a phase angle which is monotonic with respect 
to the energy” [IS] . Thus, for the Schrodinger equation, the phase angle is the logarithmic 
derivative </>(x). From equation (J2J), as x increases across a node of wave function i/)(x), 4>(x) 
decreases to —oo, jumps to +oo and then decreases again. 

In 2005, Ma and Xu [S] by carefully studying one-dimensional Schrodinger equation gen¬ 
eralized this exact quantization rule to the 3D radial Schrodinger equation with spherically 
symmetric potential by simply making the replacements x —> r and V(x) —> V e ff(r ): 

n -1 






k(r)dr = Ntt + / (j>{r) 


dk(r) 

dr 


d(j)(r) 
dr 


dr, k(r) = \1%[E -V eff (r)], (3) 


where and are two turning points determined by E = V e ff(r). The N = n + 1 is 
the number of the nodes of </>(r) in the region E n i = V e ff{r) and is larger by one than 
the number n of the nodes of wave function i^(r). The first term Nn is the contribution 
from the nodes of the logarithmic derivative of wave function, and the second is called the 
quantum correction. Ma and Xu [8] found that for all well-known exactly solvable quantum 
systems, this quantum correction is independent of the number of nodes of wave function. 
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Accordingly, it is enough to consider the ground state in calculating the quantum correction 
Qc = K(r)%dr. 

The integrals in equations (J3J) and the calculation of quantum correction term are not 
easy to obtain for various quantum mechanical problems. This motivated Serrano et al in 
2010 to propose Qiang-Dong proper quantization rule [T2] . so as to simplify the quantum 
correction terms. This rule can be summarized as follows: 

rn rrob 

/ k(r)dr = / k 0 (r)dr + mr with n = N — 1. (4) 

Jr a J roa 

In the approach, it is required to first calculate the integral on the LHS of equation (J4]) 
and then replace energy levels E n in the result by the ground state energy E 0 to obtain the 
second integral (RHS). This quantization rule has been used in many physical systems to 
obtain the exact solutions of many exactly solvable quantum systems 0 [121 US EH (22] 


3 Application to some diatomic molecular potentials 

In this section, we apply the Qiang-Dong proper quantization to study the rotation vibra¬ 
tional of some diatomic molecular potentials. Also, where necessary, we compare our results 
with the ones obtained before in the literature. 


3.1 Tietz-Wei molecular potential 


The Tietz-Wei diatomic molecular potential we examine in this section is defined as UK EH 

m 


V (r) = D 


1 _ e ~bh(r-r e ) 1 2 

1 - c h e~ bh ^- r ^ ’ 


( 5 ) 


with bh = 5(1 — Ch), r e is the molecular bond length, 5 is the Morse constant (denoted as 
fd in some other research papers), D is the potential well depth and Ch is an optimization 
parameter obtained from ab initio or Rydberg-Klein-Rees (RKR) intramolecular potentials. 
r is the internuclear distance. When the potential constant approaches zero, i.e. Ch —> 0, 
the TW potential reduces to the Morse potential ra- The shape of this potential is shown 
in Figure Q] for different molecules. To study any quantum physical model characterized 
by the diatomic molecular potential given by equation (J4]), we need to solve the following 
Schrodinger equation: 
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Figure 1: Shape of Tietz-Wei diatomic molecular potential for different diatomic molecules: (I) NO(a 4 Ih) 
(II) NO(B 2 II r ) (III) NO (L' 2 cf) (IV) NO(6 4 S-) (V) NO(V 2 n r ) (VI) H 2 (X 4 S+) (VII) CO(V 1 E+) (VIII) 

ici(v 4 e+) (ix) ici(A 3 iii) (x) ici(T' 3 n 2 ). 

i^m + ~ En £ ) ^ ^ = ( 6 ) 


In this section, we take the V(r) as the Tietz-Wei potential. Now we begin by applying 
the method of variable separation so as to split equation (EJ) into radial and angular part. 
Thus, by taking the wavefunction il} n ,e,m( r ,9, 4>) as r _1 .R n ^(r)Y£ m (#, </>) the radial part can 
be found as 


d 2 R nt (r) 2/j, 

F „ D 

1 _ g -bh(r-r a ) 

2 £(£ + l)fl 2 ' 

dr 2 h 2 


1 - C h e- b b( r ~ r e) _ 

2/ir 2 


Rni{r) = 0, 


(7) 


where n, l and E n i denote the principal quantum numbers, orbital angular momentum 
numbers and the bound state energy eigenvalues of the system under consideration (i.e., 
E n i < 0 ), respectively. It is generally known that for i = 0, problem (J7J) is exactly solvable 
but for i ^ 0, it isn’t. Therefore, in order to solve the above equation for i ^ 0 states, 
Hamzavi et al (2HJ found that the following formula 

l \ ( e ~bh(r~r e ) g-2 b h (r-r e ) 

r 2 r 2 l ° 1 1 — Che~ bh< ' r ~ re '> 2 (f _ C/i g-^h-r e )) 2 
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(a) 



c h (Dimensionless Unit) 


(b) 



c h (Dimensionless Unit) 


Figure 2: (a) The variation of the ground state energy spectrum for various values of l as a 
function of the potential constant Ch■ We choose ji = 1, bh = 5, r e = 0.8 and D = 15. (b) The 
variation of the first excited energy state for various l as a function of the potential constant Ch- 


(a) (b) 



Figure 3: (a) The variation of the ground state energy state for various values of t as a function 
of the particle mass fj,. We choose Ch = 0.03, bh = 5, r e = 0.8 and D = 15. (b) The variation of 
the first excited energy state for various t as a function of the particle mass /i. 

with 

D 0 = 1 - -(1 - c h )(3 + c h ) + -^(1 - c h ) 2 , lim D 0 = 1 — - + (9a) 

a or c h ->o a cr 

Di = -(1 - c /l ) 2 (2 + c h ) - -^(1 - c h ) 3 , lim D\ = — —(9b) 

a a z b c h ^ o a a z 

D 2 =-(1 — C/,,) 3 ( 1 + Ch) H- 2(1 — c ^) 4 ; lim D 2 = -1-2> 

a c h -i>o a a z 


(9c) 



















(a) 



(b) 



Figure 4: (a) The variation of the ground state energy state for various values of i as a function 
of the parameter bh- We choose Ch = 0.03, fi = 1, r e = 0.8 and D = 15. (b) The variation of the 
first excited energy state for various t as a function of the parameter bh- 


(a) (b) 



Figure 5: (a) The variation of the ground state energy state for various values of i as a function of 
the molecular bond length r e . We choose n = 1, bh = 5, = 0.03 and D = 15. (b) The variation 

of the first excited energy state for various £ as a function of the molecular bond length r e 


is a good approximation scheme to deal with the centrifugal potential term. Constant 
a = bhT e has been introduced for the sake of simplicity. Now, by inserting this approximation 
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(a) 


(b) 



Figure 6: (a) The variation of the ground state energy state for various values of i as a function 
of the potential well depth D. We choose g = 1, bh = 5, Ch = 0.03 and r e = 0.8. (b) The variation 
of the first excited energy state for various l as a function of the potential well depth D 


Table 1: Model parameters of the diatomic molecules studied in the present work. 


Molecules(states) 

Ch 

m/io- 23 (s) 

bRA- 1 ) 

r e (A) 

D(cm 4 ) 


Refs 

NO(a 4 n;) 

0.0082003 

1.249 

2.408413 

1.451 

16361 

2.428326 

ESI 

NO(B 2 n r ) 

-0.482743 

1.249 

3.42650 

1.428 

22722 

2.310923 

m 

NO (L ,2 </>) 

-0.073021 

1.249 

2.73796 

1.451 

14501 

2.551645 

m 

NO(b 4 S-) 

-0.085078 

1.249 

3.01538 

1.318 

21183 

2.778957 

m 

NO(x 2 n r ) 

0.013727 

1.249 

2.71559 

1.151 

53341 

2.7534 

m 

h 2 (x 1 s+) 

0.170066 

0.0837 

1.61890 

0.7416 

38268 

1.9426 

m 

CO(X 1 S+) 

0.149936 

1.1392 

2.20481 

1.1283 

9.0540 

2.2994 

m 

ICl(.Y 1 S+]) 

-0.086212 

4.55237 

2.008578 

2.3209 

17557 

1.849159 

m 

ici(A 3 ni) 

-0.167208 

4.55237 

2.542557 

2.6850 

3814.7 

2.178324 

m 

ici(A ,3 n 2 ) 

-0.157361 

4.55237 

2.373450 

2.6650 

4875 

2.050745 

m 


into equation (J7J) and then introducing a new transformation of the form r —y g = 
through the mapping function Q — fir) with r in the domain [0, oo) or g in the domain 
[—1, oo], we obtain the following second order differential equation: 


1 d 2 R ni (g) 2 n 


dg 2 


+ [Eni - Veff(e)] Rnl{g) = 0, with 
n 

£(£+l)D 0 2 gD , m $ Dl + ^(c fa * !) 


l(l+l)P2 , 2 [ID 


h 2 


+ 


+ 


+ ^(c h ~ l ) 5 


(e“» - c h y 


( 10 ) 

2/i 


K//(e) 


e<*e - c h 


























Table 2: Comparison of the bound-state energy eigenvalues —(E n i — D){eV) of th and CO 
molecules for various n and rotational l quantum numbers in Tietz-Wei potential. 


n 

l 

PQR (Present) 

h 2 

NU (Present) 

NU [24] 

GPS [30] 

PQR (Present) 

CO 

NU (Present) 

NUl24l 

0 

0 

4.48160205 

4.48157183 

4.4815718267 

4.4815797825 

11.07378271 

11.07370964 

11.07370964 


5 

4.26706288 

4.26703273 

4.2658220403 


11.06668238 

11.06660933 

11.06659606 


10 

3.73741847 

3.73738873 

3.7336304360 


11.04775337 

11.04768033 

11.04763173 

5 

0 

2.28155518 

2.28153387 

2.2669650930 

2.2815913849 

9.63231601 

9.63224796 

9.629868985 


5 

2.12184477 

2.12182368 

2.1070072990 


9.62560750 

9.62553946 

9.623148913 


10 

1.72567170 

1.72565129 

1.7105026080 


9.60772327 

9.60765524 

9.084920084 

7 

0 

1.63001750 

1.62999964 

1.6130911000 


9.08816748 

9.08810144 

9.084920084 


5 

1.48916801 

1.48915041 

1.4722799590 


9.08161259 

9.08154656 

9.078354525 


10 

1.13910018 

1.13908334 

1.1225356530 


9.06413799 

9.06407196 

9.060851440 


Table 3: Bound-state energy eigenvalues for NO (a 4 IIi), NO (T? 2 n r ), NO (L ,2 </>), NO (6 4 E _ ) 
and NO (A" 2 n r ) molecules for various n and rotational £ quantum numbers in Tietz-Wei 
diatomic molecular potential. 


n 

t 

NO (a 4 n ; ) 

NO (B 2 n r ) 

NO (L' 2 0) 

NO (fc 4 E-) 

NO (x 2 n r ) 

0 

0 

0.05724382 

-0.06511639 

-0.05748847 

-0.07561426 

-0.04775936 

1 

0 

0.16924769 

-0.19664818 

-0.17491875 

-0.22970681 

-0.14364066 


1 

0.16950472 

-0.19637388 

-0.17464813 

-0.22937970 

-0.14352509 

2 

0 

0.27795077 

-0.33036890 

-0.29569047 

-0.38771825 

-0.24006757 


1 

0.27820318 

-0.33009351 

-0.29541545 

-0.38738640 

-0.23995169 


2 

0.27870799 

-0.32954274 

-0.29486542 

-0.38672273 

-0.23971993 

3 

0 

0.38335537 

-0.46625489 

-0.41978250 

-0.54962300 

-0.33703799 


1 

0.38360316 

-0.46597847 

-0.41950311 

-0.54928644 

-0.33692180 


2 

0.38409872 

-0.46542563 

-0.41894434 

-0.54861334 

-0.33668942 


3 

0.38484202 

-0.46459641 

-0.41810624 

-0.54760374 

-0.33634086 

4 

0 

0.48546381 

-0.60428372 

-0.54717398 

-0.71539576 

-0.43454981 


1 

0.48570698 

-0.60400632 

-0.54689024 

-0.71505453 

-0.43443332 


2 

0.48619332 

-0.60345153 

-0.54632279 

-0.71437207 

-0.43420033 


3 

0.48692277 

-0.60261938 

-0.54547165 

-0.71334844 

-0.43385085 


4 

0.48789530 

-0.60150990 

-0.54433688 

-0.71198366 

-0.43338490 


The two turning points are obtained by solving V e ff(g) — E n e = 0 or V e ff(p) — E n i = 0 with 
p = (e ag — Ch ) -1 . Thus, it is easy to show that the turning points p a and pb are 


‘^ + ^(c h -l) + J\4^ + mc h ~l) 


Pa = 


— AT, 


nt 


+%?(<* - 1 )* 


2 T, 


nl 


(11a) 
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Table 4: Bound-state energy eigenvalues for H 2 (X 1 !!^), IC1 IC1 (k^fli), IC1 

(v 4 /3 II2) and IC1 ( B'0 + ) molecules for various n and rotational i quantum numbers in 
Tietz-Wei diatomic molecular potential. 


n 

l 

h 2 (x 1 e+) 

CO (X 1 E+) 

IC1 (X 1 E+) 

IC1 (/Fill) 

IC1 (A' 3 n 2 ) 

0 

0 

0.26925518 

0.15070333 

-0.02388751 

-0.01317364 

-0.01400330 

1 

0 

0.77877437 

0.44839050 

-0.07201091 

-0.03995431 

-0.04239940 


1 

0.79255794 

0.44885861 

-0.07198238 

-0.03993286 

-0.04237766 

2 

0 

1.25241802 

0.74134891 

-0.12061278 

-0.06734751 

-0.07134441 


1 

1.26542046 

0.74181176 

-0.12058411 

-0.06732587 

-0.07132249 


2 

1.29131118 

0.74273745 

-0.12052676 

-0.06728259 

-0.07127867 

3 

0 

1.69133322 

1.02961057 

-0.16969191 

-0.09534658 

-0.10083355 


1 

1.70358901 

1.03006820 

-0.16966309 

-0.09532475 

-0.10081146 


2 

1.72799429 

1.03098343 

-0.16960544 

-0.09528108 

-0.10076729 


3 

1.76433953 

1.03235625 

-0.16951898 

-0.09521559 

-0.10070103 

4 

0 

2.09661354 

1.31320715 

-0.21924707 

-0.12394503 

-0.13086213 


1 

2.10815488 

1.31365958 

-0.21921810 

-0.12392301 

-0.13083988 


2 

2.13113862 

1.31456443 

-0.21916017 

-0.12387897 

-0.13079536 


3 

2.16536983 

1.31592167 

-0.21907327 

-0.12381290 

-0.13072858 


4 

2.21056327 

1.31773125 

-0.21895740 

-0.12372482 

-0.13063955 


+ ¥(<=» -1) 


^ -1) 


1 2 


Pb = - 

with the following sum and product properties: 

^ + l)D l h 2 + Aiir 2 e D{c h -l) 


4 T, 


ni 


•msi + ^(c b - i )^ 


2 T, 


ni 


Pa T Pb 


and p a p h = 


(lib) 

’+l)D 2 h 2 + 2pr 2 e D(c h -l) 2 


T n ,rl% 2 T ni r* e h 2 

Furthermore, the momentum k(p) between two turning points can be found as: 


( 12 ) 


Kp) = 


£(£+l)D 2 2 pD 

--a-+ -rr( c h - l ) 2 


[(pb - p)(p - pa)}- 


(13) 


The Riccati relation given by equation (J2J) can be re-written for the ground state as 


apil+CHp) ^ + Mp)2 = _ 2 ^ [£ , m _ K/A/>)] 
' e ft 


(14) 


Since the logarithmic derivative 0o (p) for the ground state has one zero and no pole, it has 
to take the linear form in p. The only possible solution satisfying equation (TT4l) is of the 
form 0o (p) = A + Bp. The substitution of this expression into equation (JHJ), one has the 
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ground state energy eigenvalue 


Eqi — 


£(l+l)D 0 A , 


rt 


-h D with A = —— 

2/i 2£ 


A£ + i> Dl + <£(*_!) 


r 


+ 


a 


and B = ^-^,/l + 4r « 


2ry 2rv- 


2^2 


a'c 


h 


ArE + -1) 2 

r? ft 2 


2r e 

(15) 


We have now reached a position of calculating the integrals given by equation ([!]) . The LHS 
integral can be calculated as follows: 

rPa k(p) 


rr B 


n 8 b 


k(r)dr — r e I k(g)dg = —1 


'r A 


' 8a 


« J Pb Pi 1 + c hP) 


dp 


Te 

a 


+ 1)D 2 2 pD 


+ _ ^2 _ ( Ch ~ 1)' 


rPb V(Pb~ P)(P~ Pa) 


'Pa 


P(1 + c hP ) 


dp 


(16 


7rr e £(£ + l)D 2 2pD 


a 




\/(l + c h pa){ 1 + c h p b ) 1 


Ch 


Ch 


\/7>aPb 


nr P 


a 


Trtf - 


Ch 


+ l)-^i , ^PE ^ lN 

o -t 2 ( C /i 1) 




\fRx 


+ ~ 


V^nl 


with T n£ = 


+ 1 )- p o , ?E (D _ F \ 
„2 + -*-2 \E E n i) 


tf 


and Rt = 


C h C h 

£{£ + l)D 2 h 2 


2 fir 


+ D(ch — 1)" 


where we have utilized the properties (I12p and the integral relation given by 
f Xb y/(x b -x)(x- x a ) 


x(l + Qx) 


-dx = 7r 


a / (Qx a + 1 )(QXb + 1 ) 1 , - 

- Q - Q-V^i 


(17) 


Now, simply by replacing E n £ in the above equation (TT6]i by E 0 i given by equation (lT5]h and 


T n i as Tqi 


^ + p D ° + |y (D — E 0 i) , we obtain the integral in the RHS of equation (jTJ) as 


rrob 


k 0 (r)dr = 


7 T 7 V 


"r Oa 


a 

nr e 


Toe - 


Ch 


rf fi 


Rt 
H—2“ 

Ch 


a/ Rt 


Ch 


OLCh 


\J Rt + B 


\/Toi 

(18) 


With equations (THU) . (ITS]) and (j3J), we can deduce the following relation 


irr P 


a 


T n £ - 


Ch 


±iihi + <£(* -1) 

rj h 


, ^ /tT" 

H-2- 


-B = nir. 


OLCh 


T n i - 


Ch 


+ T)D x 4 pD 

~3 - ~T2~\ Ch ~ 1) 


n z 


Ri 


+ - t - = -( n +—‘ 8 )- + Vt ^ ( 19 ) 

' ac h J r e 
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On squaring up both sides of equation (fT9|) . it is straightforward to show that the energy- 
eigenvalues equation can be found as 


Eni — 


2*2 


2 nr, 


1)A) ! D 


2 //rl 


rf + ^(PlO, - g 2 ) + ^ (l - ^ 
2rj 


( 20 ) 


1 1 


with v = n+ - + -\jl + — 


4 (D 2 £(£ + 1) 2 iiDr\ 


a 


+ ^( i -oq. 


3.2 Shifted Deng-Fan molecular potential 


The shifted Deng-Fan molecular potential model we examine in this section is defined as 

HZlEZj 

V(r) = d(i--^—\ -D, b = <?"■- 1, (21) 

where (D,D), b and /3 are three parameters representing the dissociation energy, the po¬ 
sition of the minimum r e and the range of the potential respectively. Very recently, Wang 
and co-workers found that the Manning-Rosen, Deng-Fan and Schioberg potential are not 
better than the traditional Morse potential in simulating the atomic interaction for diatomic 
molecules [28] - In order to overcome this problem, Hamzavi et ah suggested a modification 
to the Deng-Fan potential, which they referred to as the shifted Deng-Fan potential (sDF) 
mi. This modification is simply a Deng-Fan potential [191120] shifted by dissociation energy 
D ■ The researchers [iTT] examined the Schrodinger equation with this potential and ap¬ 
plied their results to some diatomic molecules na From their plot for the shifted Deng-Fan 
potential and the Morse potential using the parameters set for H 2 diatomic molecule, it was 
shown that the two potentials are very close to each other for large values of r in the regions 
r « r e and r > r e , but they are very different at r 0. Also, if both the Deng-Fan and the 
shifted Deng-Fan potentials are deep (that is, D » 1) they could be well approximated by 
a harmonic oscillator in the region r « r e dt 

In Figure (JZJ), we study the variation of this potential with respect to some diatomic 
molecules of interest given in table [Q Now inserting this potential into the Schrodinger 
equation, and then use the approximation of the form m 


,-pr 


( 1 -e-^rJ 


( 22 ) 
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Figure 7: Shape of shifted Deng-Fan diatomic molecular potential for different diatomic molecules: 
(I) NO(a 4 IF) (II) NO(B 2 U r ) (III) NO(L ,2 </>) (IV) NO(6 4 £-) (V) NO{X 2 H r ) (VI) H 2 (V 4 £+) (VII) 

co(v 1 e+) (viii) ici(v 1 e+) (ix) ici(^ 3 ni) (x) ici(T , 3 n 2 ). 


(a) (b) 



Figure 8: (a) The variation of the ground state energy state for various values of t as a function 
of the potential range /?. We choose fi = 1. r e = 0.8 and D = 15. (b) The variation of the first 
excited energy state for various £ as a function of the potential range /3. 
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(b) 



U (a.u.) 



Figure 9: (a) The variation of the ground state energy state for various values of £ as a function 
of the particle mass y. We choose /3 = 5, r e = 0.8 and D = 15. (b) The variation of the first 
excited energy state for various £ as a function of the particle mass y. 


(a) (b) 



Figure 10: (a) The variation of the ground state energy state for various values of £ as a function 
of the molecular bond length r e . We choose y = 1, /3 = 5 and D = 15. (b) The variation of the 
first excited energy state for various £ as a function of the molecular bond length r e 


the effective potential takes the following form: 

V e ff(y) = P + Qy + Ry 2 with P = D - D + 


£{i + 1 )(3 2 d 0 h 2 

2/i 


(23) 


_ t(t + \ )/? 2 a 2 , _ ne + i)p 2 d„h 2 2 

Q = — -2 Db anf| R = — ^ h Db 2 , 


2/i 


2/i 
























(a) (b) 



Figure 11: (a) The variation of the ground state energy state for various values of £ as a function 
of the dissociation energy D e . We choose // = 1, /3 = 5 and r e = 0.8. (b) The variation of the first 
excited energy state for various £ as a function of the dissociation energy D 


Table 5: Comparison of the bound-state energy eigenvalues —E n i(eV) of H 2 and CO molecules 
for various n and rotational £ quantum numbers in sDF diatomic molecular potential. 


n 

1 

H 2 Diatomic Molecule 

Present GPS [2g] AIM [27] 

n-u |T7] 

Present 

CO Diatomic Molecule 

GPS Hi] AIM [27] N-U[23| (CO) 

0 

0 

4.39461978 

4.39462330967 

4.394619779 

4.39444 

11.08074990 

11.0807513815 11.08075178 

11.08068 


5 

4.17661316 


4.176618048 

4.17644 

11.07253746 

11.07253985 

11.07247 


10 

3.62182049 


3.621838424 

3.62165 

11.05064208 

11.05064581 

11.05057 

5 

0 

1.75845157 


1.758451567 

1.75835 

9.68814442 

9.688146187 

9.68809 


5 

1.61740572 


1.617410615 

1.61731 

9.68022402 

9.680226284 

9.68017 


10 

1.26043371 


1.260451640 

1.26034 

9.65910731 

9.659110919 

9.65905 

7 

0 

1.07763699 


1.077636993 

1.07756 

9.15916229 

9.159164003 

9.15911 


5 

0.96180989 


0.961814782 

0.96174 

9.15135744 

9.151359661 

9.15131 


10 

0.66982613 


0.669844065 

0.66976 

9.13054886 

9.130552425 

9.13050 


after an appropriate coordinate transformation of the form y = [eP r — l) 1 has been intro¬ 
duced. Now, we can write the non-linear Riccati equation for the ground state as 

-ay{ 1 + yW Q {y) + <j>l(y) = -p [Ew - V eff (y)\ (24) 

Since the logarithmic derivative <j>o(y) for the ground state has one zero and no pole, it has 
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Table 6: Bound-state energy eigenvalues for NO (a 4 !!*), NO ( B 2 U r ), NO (L ,2 </>), NO (6 4 £“) 
and NO (X 2 II r ) molecules for various n and rotational £ quantum numbers in sDF diatomic 
molecular potential. 


n 

i 

no (o 4 ni) 

NO (b 2 n r ) 

NO {l' 2 <£) 

NO (fe 4 xr) 

NO (v 2 n r ) 

0 

0 

-1.96954814 

-2.75045354 

-1.73990613 

-2.54988834 

-6.49084320 

1 

0 

-1.85431702 

-2.61950265 

-1.62684906 

-2.40042016 

-6.24919850 


1 

-1.85394591 

-2.61914169 

-1.62645839 

-2.39995170 

-6.24866000 

2 

0 

-1.74265661 

-2.49186027 

-1.51767759 

-2.25557591 

-6.01232867 


1 

-1.74228862 

-2.49150219 

-1.51729002 

-2.25511093 

-6.01179390 


2 

-1.74155263 

-2.49078605 

-1.51651490 

-2.25418099 

-6.01072438 

3 

0 

-1.63455757 

-2.36751730 

-1.41238208 

-2.11534494 

-5.78022207 


1 

-1.63419268 

-2.36716213 

-1.41199768 

-2.11488343 

-5.77969103 


2 

-1.63346278 

-2.36645179 

-1.41122876 

-2.11396043 

-5.77862896 


3 

-1.63236803 

-2.36538629 

-1.41007542 

-2.11257593 

-5.77703588 

4 

0 

-1.53001045 

-2.24646485 

-1.31095323 

-1.97971663 

-5.55286713 


1 

-1.52964860 

-2.24611256 

-1.31057185 

-1.97925858 

-5.55233981 


2 

-1.52892485 

-2.24540801 

-1.30980918 

-1.97834250 

-5.55128517 


3 

-1.52783925 

-2.24435119 

-1.30866514 

-1.97696839 

-5.54970324 


4 

-1.52639185 

-2.24294214 

-1.30713977 

-1.97513628 

-5.54759404 


to take the linear form in y. Thus, we assume the following solution for the ground state 


Mv) = A + By. 


(25) 


By putting equation (1251) into (1351) and then solve the non-linear Riccati equation, it is 
straightforward to obtain the ground state energy and values of A and B as 


Em — P — 


h 2 A 2 

2y 


with A = 


/i Q 


nr 


- R B 
~B + 2" 


and 


B = ^ + i 
2 2 


P 2 + 


8yR 


(26) 


Furthermore, in a similar fashion to the previous problem, the two turning points as well as 
their sum and product properties are given by 


Va = ~§R~lR^ Q2 ~ 4R(P ~ Ent) ’ ™ d Vb= ~§R + ~ 4R{P ~ En,) 


Va + Vb ~ 


Q 

R' 


P - Pni A ; / \ 
VaVb = - R - and k ^ = 


2 fiR 


[- (y - Va) (y - w)] 1/2 ■ (27) 


Now, we have all necessary tools required to perform our calculations. Therefore, we proceed 
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Table 7: Bound-state energy eigenvalues for H 2 (X 1 !!^), IC1 (X 1 !!^), IC1 (A 3 IIi), IC1 
(v4 /3 ri2) and IC1 ( B'0 + ) molecules for various n and rotational ^ quantum numbers in sDF 
diatomic molecular potential. 


n 

l 

h 2 ( x 3 e +) 

CO (X 1 S+) 

ici ( x 3 e +) 

ICI (u 3 ni) 

ici (u ,3 n 2 ) 

0 

0 

-4.40110165 

-11.08006149 

-2.152719356 

-0.459938084 

-0.59052741 

1 

0 

-3.75459356 

-10.79430704 

-2.104958026 

-0.434430900 

-0.56321617 


1 

-3.74109390 

-10.79376640 

-2.104907368 

-0.434368596 

-0.56316009 

2 

0 

-3.17123615 

-10.51257523 

-2.057739707 

-0.409653282 

-0.53655412 


1 

-3.15882350 

-10.51203847 

-2.057689131 

-0.409591035 

-0.53649810 


2 

-3.13406357 

-10.51096496 

-2.057587980 

-0.409466541 

-0.53638607 

3 

0 

-2.64736026 

-10.23485368 

-2.011064148 

-0.385605052 

-0.51054106 


1 

-2.63597228 

-10.23432079 

-2.011013655 

-0.385542864 

-0.51048511 


2 

-2.61325755 

-10.23325500 

-2.010912669 

-0.385418487 

-0.51037321 


3 

-2.57933793 

-10.23165635 

-2.010761189 

-0.385231921 

-0.51020535 

4 

0 

-2.17959586 

-9.96113006 

-1.964931100 

-0.362286035 

-0.48517681 


1 

-2.16917525 

-9.96060101 

-1.964880689 

-0.362223905 

-0.48512092 


2 

-2.14839142 

-9.95954293 

-1.964779867 

-0.362099645 

-0.48500914 


3 

-2.11735861 

-9.95795584 

-1.964628634 

-0.361913255 

-0.48484148 


4 

-2.07624689 

-9.95583977 

-1.964426991 

-0.361664734 

-0.48461792 


to calculate integral (j3J) 



cyb 


Hv) 


!y a Py{ 1 + y) dy Jy a 


Vb 2 yR[(y-y a )(y b -y)} 1/2 


it 12yR 

1 


71 /2{lR 

T 


/ 3 2 h 2 y(i + y) 
\/(i + y a )(i + yb) - l - 
y/R — Q + P — E n i 


dy 


(28) 


R 


P — E, 
R 


'nl 


where we have used the following standard integral 


ry b 


'Va 


(V ~ 2/a) (y ~ Vb)} 

2/(1 + y) 


1/2 


= 71 


Vida + !)(//&+ 1) - 1 - y/Vdjb 


(29) 
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Furthermore, we can find 


rr ob 


ko(r)dr = 


fVOb 


k(y) 


'roa 


'VOa 


PvO- + y) 


dy = 


fVOb 


'VOa 


2 yR [- (y-y a ) {y-yb)] 1/2 


a 2 h 2 


2/(1 + y) 


71 2 flR 


p\ h 


V (! + 2/a)(l + Vb) - 1 - VVc 


yb 


(30) 


7r 2 yR 


/9 V h 2 


y/R ~ Q + P — Eoe _ ^ _ IP ~ ^ ol 


n 2 fiR 


f3\ h 2 


h 2 


2/iR 


R 


(■ A-B ) 


R 


P - E oe 

71 

/2/j/? 

B\ 


R 

~P\ 

* h 2 

\ 

2/i,/? 


From equations (j3J), (I29|) and Oil) , we can find the energy spectrum for the sDF as 


2/i 


2p 


2h j3 


ZR2 I b 


+ n 


- D. 


(31) 


4 Calculation of the eigenfunctions 


Eigenfunctions-eigenvalue relation is very important in quantum mechanics because of its 
prominence in the equations which relate the mathematical formalism of the theory with 
physical results, eigenfunctions could be considered as trial functions in variational-type 
procedures for deriving energy levels anl also for computing line intensities. Since proper 
quantization rule cannot be used to obtain these eigenfunctions, we therefore resort to using 
the recently proposed formula method (Tj. This method is very easy to use in obtaining 
not only the eigenfunctions but also energy eigenvalues. In the approach, it is required 


to transform the Schrodinger equation with two solvable quantum molecular systems-Tietz- 
Wei and shifted Deng-Fan potential models into the form given by equation (1) of ref. [lj via 
an appropriate coordinate transformation of the form r = e^ e (for TH) and t = e~ ar (for 
sDF), which maintained the finiteness of the transformed wave functions on the boundary 


conditions to have 

d 2 R n£ (T) 1 dR nl (r) 


dr 2 
-2 c h 


t dr 
2 lir 2 e E, 


ni 


+ 


a 2 h 2 
2 yr 2 e E ni _ 
a 2 h 2 


r 2 (l - C h T) 2 

+ 1 ) 




-Di 


cr 

+ 1) 


h 2 a 2 
Anr 2 D 


cr 


+ 1 ) 


Dr 


CT 


h 2 a 2 

+ 1 ) 


cr 


D 1 


T 


CT 


m n ^ 2 nr 2 e D 

{DiC h - D 2 ) 2 2 

n a z 
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Rni(r) = 0 , 

(32a) 


























































d 2 R n e{t) 1 dRjfy 1 p fi ‘liiDb (A fibD 


—2 (En, — D) —1(1 + 1)(1 - 2do)') + « 2 f Ap (£„, — D) —1(1 + 1H 


fln<(t) = 0. 
(32b) 


Considering equation (I32al) with reference to P, k\, k 2l k%, At-h, Bt-h and Ct-h can be 
found. Then, parameters k 4 and k$ can be obtained as 


ki = 


^ (D _ Eni) + ^ Dn 


h 2 a 2 


a “ 


, , 1 I / 4 

and k 5 = -\ l + \ll + ~2 


h - 


'+ 1 ) n , 2 , 9 
- B 2 + j. 9 „ (1 — Cft) 2 


cr 


h 2 a 2 


(33) 


Hence, the eigenfunctions for TW can be found as 


Rnt(o) = N ni e k ± ae (l - Ch e ae ) ks 2 Fi (-ra, n + 2 (k 4 + k 5 ); 2 k A + 1, c h e ae ) . (34) 

Similarly, the eigenfunctions for sDF can be found as 


Rni(z) = N ne t w (l - t) v 2 Fi (-n, 2 (w + v) -,2w + l] t ), 


(35) 


where 


w = 


f (£W - B) - f^(!» + 2) + 


+ l)d„) and t,-i + pC+ip^-. 

(36) 


5 Numerical Results and Discussion 

In FigureP we plotted the Tictz-Wei (TW) potential for different diatomic molecules. In 
what follows, to see the behavior of the ground n = 0 and first excited n = 1 states, we 
plotted the energy for these states with potential parameters for three different orbital states 
i = 0,1,2. In Figure El we show the variation of E n ^ with the potential constant Ch■ It 
shows that for Ch < 0, the energy is negative whereas when Ch > 0, the energy is positive 
for n = 0. On the other hand, for n — 1, the energy becomes strongly bound for Ch < 0 and 
moves toward the negative energy for Ch > 0. The Ch = 0 represents the Morse energy. The 
best choice C/,, = 0.03 restores the results of Morse potential. At this value the energy curves 
coincide and have same behavior for £ = 0,1,2. Figure [3] shows the variation of E n ^ with 
the reduced mass /1 for three orbital states. The energy is very similar for 0.1 < /i < 1.0 
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when n — 1 but different when n = 0. Its seen that when /i increases for more than 0.3, the 
energy spectrum becoming positive for ground state while for excited state, it is positive for 
any value of /i 

In Figure 0 we plotted the variation of E n< g with the potential parameter bh- It is 
increasing in the positive direction within the interval 0 < bh < 8 when n — 0. However, 
when n — 1, the energy increases in positive side for 0 < bh < 4 and increases in the negative 
side for 4 < b h < 8. In FigureOwe show the variation of the energy states E n ^ as a function 
of molecular bond length r e . The groimd state energy drops with nearly 2eV for all orbital 
states at r e = 0.8 fm and r e = 0.9/m whereas the first excited state has drop of about 
0.45eV and coincide at 0.56fm. Finally, Figure El demonstrates the energy versus the well 
depth D. Its seen that the ground state energy span from negative to positive spectrum at 
D = 2eV. for orbital states i = 0,1,2. However the first excited energy state span from 
negative to positive spectrum at D — 8eV for all orbital states. 

A very similar behavior to TW potential model (sDF shape) for various molecules 
is shown in Figure 7. In addition we have obtained the energy spectrum for different 
diatomic molecules with the help of TW molecular model for various states using the 
model potential parameters in Table |T] This spectroscopic parameter are taken from Refs. 
m, ra and [31] and the conversion factors used are taken from NIST database [32]: 
1 cm ” 1 = 1.239841930eV, he = 1973.29e1/A and 1 amu = 931.494061Meu/c 2 . In Table [2] 
we test the accuracy of the method utilized in this study by finding the energy spectra of H 2 
and CO diatomic molecules. We found that the spectrum obtained by NU method in [23] 
have some error in the Maple codes. We therefore re-compute these spectrum in the present 
work for the sake of comparison. As it can been seen from the table, our results are very 
close to the ones of the Nikiforov-Uvarov method. Tables [3] and [4] present the spectrum for 
H 2 CO (X 1 E + ) and various electronic states of NO and IC1 diatomic molecules. 

Considering sDF molecular potential, Figure [S] shows the variation of E n ( i as a function 
of (3, in the ground state. The restriction on choice of the parameter /3 of sDF molecular 
potential can be observed. The energy for i = 1,2 increases in the positive side but i = 0 
the energy increases in the negative side. On the other hand, in the first excited state, the 
energy increases in the negative side in the interval 0 < /3 < 15 for i = 0,1, 2. 

Figure [9] shows the variation of E n ^ as a function of reduced mass p. Small values of 


20 


particle mass /i result into a sharp change in energy values for ground and first excited states 
for the orbital states. The energy becomes stable when // > 1. This plot indicates how to 
choose or read the most reasonable values of /i which provides the most appropriate and 
not overlapping spectrum amongst orbital states. Figure [TO] is a plot of energy versus bond 
length r e . The spectrum > E 0t 1 > E 0)0 when r e > 0.2 and E\- 2 > E\j > Ei t0 when 
r e > 0.4. Figure 10 set restrictions on the most suitable values of r e . At r e < 0.6 — 0.8/m, 
the energy of different orbital states overlap and deteriorate sharply. 

The variation of the energy versus parameter D is shown in Figure [TTJ For n = 0 the 
energy increases and then decreases in the given range 0 < D < 25 whereas for n = 1, it is 
increasing in the same interval. Furthermore we generated the spectrum of several diatomic 
molecules using the sDF molecular potential for various states. The behavior of the plot 
energy against each potential parameter for various states provides us the most appropriate 
physical values for each parameter. Figures 6 and 10 have different behaviors since they 
are a plots of energy against r e for two potentials. However, if one has set to choose large 
values for r e (say r e > 0.6) then the two curves will be similar. 

Table O compares our results for H 2 and CO with those of the GPS method, Nikiforov- 
Uvarov method and AIM methods. Our currently found energy states are reasonably com¬ 
pared with the other findings. The vibrational energy is close to 7 digits with AIM [27] 
but found to agree with GPS [26] up to 4 digits. However, the rotational-vibrational en¬ 
ergy states are close close to 5 digits with AIM. This is due to the approximation made 
to the centrifugal restorsion term. Also, it should be noted that the model is a parameter 
dependent which may result into slight variation in energy spectrum if parameters are not 
converted/adjusted properly. Table [6] displays the energy spectrum for different species 
of NO diatomic molecules. However, Table [7] present ones for H 2 CO (2T 1 £ + ), 

various electronic states of I Cl diatomic molecules 


6 Concluding Remarks 

In this research work, we applied proper quantization rule in a spectroscopic study of some 
diatomic molecules. This task is made possible by solving the Schrodinger equation with two 
molecular models; namely, Tietz-Wei and shifted Deng-Fan potential models. This solution 
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serves as the basis for the description of the quantum aspects of diatomic molecules. We 
obtained the energy spectra of different diatomic molecules. The validity and accuracy 
of the method is tested with previous techniques via numerical computation for H 2 and 
CO molecules. Our reasonable results show the efficiency and simplicity of the present 
calculations. The approximation to the centrifugal restorsion is valid for the lowest orbital 
quantum number l. As £ increases, the accuracy of the energy states reduces and vice- 
versa. The present research work represents a new procedure in dealing with the diatomic 
molecules. Our results are reasonable and credible in generating the spectrum as the other 
commonly known methods. 
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